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Abstract. We present results on the chiral and deconhnement properties of the 

QCD transition at finite temperature. We performed calculations using asqtad and 

HISQ/tree actions using lattices with temporal extent N T = 6, 8 and 12 allowing to 

£NJ \ control the approach to the continuum limit. We analyze the chiral transition in terms 

of universal O(N) scaling functions. From peaks in the scaling functions we perform a 
simultaneous continuum extrapolation for HISQ/tree and asqtad to derive the critical 
temperature, T c = 157 ± 6 MeV. 



1. Introduction 



f^. | Improved staggered fermion actions are widely used to study QCD at non-zero 



temperatures and densities, see e.g. Ref. pQ. However, discretization effects for 



(N 

O 

lO \ staggered fermion actions are quite large in the low temperature region. To control the 

discretization effects it is important to perform calculations for large values of temporal 

extent N T or use actions with smallest possible discretization effects. Here we report 

the study of the chiral and deconfinement aspects of the finite temperature transition 

using the asqtad and HISQ/tree action on lattices with temporal extent N T = 6, 8 and 

12 with light quark masses mi = m s /20 with m s the physical strange quark mass. This 

corresponds to the lightest pion mass of about 160 MeV. The lattice spacing was set 

using the static potential. The details of the lattice setup and analysis were in part 

presented in j2] and will be discussed in a forthcoming publication [3] . 

2. Chiral transition 

For vanishing light quark masses there is a chiral phase transition which is expected to 
be second order and thus governed by universal 0(4) scaling. However, even for non- 
vanishing light quark masses, provided they are small enough, universal scaling allows 
to define pseudo-critical temperatures for the chiral transition. Thus when studying the 
chiral transition one first needs to establish that the lattice results can be described 
in terms of 0(4) scaling. In the staggered fermion formulation there is one further 
complication. Since this formulation only preserves a part of the chiral symmetry, the 
relevant universality class in the chiral limit for non-vanishing lattice spacing is actually 



Study of chiral and deconfinement transition in lattice QCD with improved staggered actions2 




HISQ/tree: N T = 8 

m/m s = 0.050, 0(4) 

0(2) 

0.025,0(4) 

0(2) 

0.037,0(4) 

0(2) 



HISQ/tree: N T = 12 



rri|/m s = 0.050, 0(4) 




140 150 160 170 180 190 200 210 220 



Figure 1. The numerical results of the order parameter as function of the temperature 
obtained with HISQ/tree action on N T = 8 (left) and N T — 12 (right) lattices together 
with the scaling fits. 



0(2). Fortunately, in the numerical analysis the differences between 0(2) and 0(4) 
universality classes are small so when referring to scaling we will use the term O(N) 
scaling. Previous studies with the p4 action provided evidence for O(N) scaling [HIT]. 
Therefore we would like to check if O(N) scaling works for asqtad and HISQ/tree action. 
In Fig. [T]we show our numerical results obtained with HISQ/tree action for the order 
parameter 



M h = 



rp4 



(1) 



For N T = 6 and N T = 8 we also made use of preliminary results at smaller quark mass 
mi = m s /40 [4]. We performed a fit of M b to the scaling form 

M b (T, mi, m.) = h 1/5 f G (t/h 1/ps ) + f M ,re g (T, H). (2) 

Here we defined H = mi/m s and the reduced temperature t and the magnetic field h 




h = -H, 

h 



(3) 



T c ° is the critical temperature in the chiral limit. Furthermore, fM,reg{T, H) = (a t AT + 
b\)H parametrizes the contribution from the regular part of the free energy density. By 
performing a 5 parameter fit of the numerical results on M b , i.e. treating T^,t ,h ,a t 
and b\ as fit parameters we can describe the temperature and quark mass dependence 
of Mb. Both 0(4) and 0(2) scaling fits work well. The pseudo-critical temperature can 
be defined as peak position of the chiral susceptibility 

d 



Xm,l 



dmi 



W)i , q = l, s . 



(4) 



Since the scaling Ansatz describes the quark mass dependence and temperature 
dependence of M b and Xm,i it can be used to determine the peak positions for different 
N T for the physical value of the light quark mass mi/m s = 27.3. Then performing a 
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combined 1/N% extrapolation of T c values obtained with asqtad and HISQ/tree action 
as shown in Fig. 2 we obtain 

T c = (157±4±3±1) MeV, (5) 

where the first error is statistical, the second error is the systematic error due to 
continuum extrapolations and uncertainties related to the differences in 0(2) and 0(4) 
scaling fits, and the last error is the overall error on the lattice scale determination. To 
present a combined error we add the first two errors in quadrature and then add the 
third one. This gives T c = 157 ± 6 MeV. 

3. Deconfinement transition 

The Polyakov loop is an order parameter for the deconfinement transition in pure gauge 
theory, which is governed by Z(N) symmetry. For QCD this symmetry is explicitly 
broken by dynamical quarks. There is no obvious reason for the Polyakov loop to be 
sensitive to the singular behavior close to the chiral limit although speculations along 
these lines have been made [5]. The Polyakov loop is related to the screening properties 
of the medium and thus to deconfinement. After proper renormalization, the square 
of the Polyakov loop characterizes the long distance behavior of the static quark anti- 
quark free energy; it gives the excess in free energy needed to screen two well-separated 
color charges. The renormalized Polyakov loop has been studied in the past in pure 
gauge theory [HI [9] as well as in QCD with two [TO], three [H] and two plus one flavors 
[I~2~| H3"] . The renormalized Polyakov loop, calculated on lattices with temporal extent 
N T , is obtained from the bare Polyakov 

/l Nt ^ \ 

L ren (T) = z(f3) N ^L bare (f3) = z(f3) N ^ ( -Tr J] U (x ,x) ) , (6) 

\ 6 x =0 I 

where z((3) = exp(— c(/3)/2) and c(/3) is the additive normalization of the static potential 
chosen such that it coincides with the string potential at distance r = 1.5r with r 
being the Sommer scale. The numerical results for the renormalized Polyakov loop for 
the HISQ/tree action are shown in the right panel of Fig. [2] as function of T/T c . As 
one can see from the figure the cutoff (N T ) dependence of the renormalized Polyakov 
loop is small. We also compare our results with the continuum extrapolated stout 
results [15] and the corresponding results in pure gauge theory [8j |9]. We find good 
agreement between our results and the stout results. We also see that in the vicinity of 
the transition temperature the behavior of the renormalized Polyakov loop in QCD and 
in the pure gauge theory is quite different. 

4. Conclusions 

We have studied the chiral and deconfinement aspects of the finite temperature 
transition in QCD. The chiral pseudo-critical temperature defined as peak of the chiral 
susceptibility for the physical quark masses was found to be 157 ± 6 MeV. The chiral 
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Figure 2. The continuum extrapolation of T c (left) and the renormalized Polyakov 
loop as function of T/T c (right). For stout data we used the value of T c — 157 MeV 
from the inflection point of the renormalized chiral condensate [15] . 



transition temperature obtained in Refs. [14] 115] using different observables was found 
to be in the range 147 — 157 MeV in good agreement with our result. Previous attempt 
to determine the chiral transition temperature by RBC-Bielefeld collaboration resulted 
in too high value of T c = (192±7±4) MeV because no reliable continuum extrapolation 
can be done by using only N T = 4 and 6 lattices. The renormalized Polyakov loop does 
not show rapid change in the vicinity of the chiral transition temperature and is quite 
different from the pure gauge theory result. Our present findings are very similar to the 
previous results obtained using the stout action [HI [15] . 
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